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Abstract 

We study the response of a single degree of freedom mechanical system com- 
posed of a primary mass, M, a linear spring, a viscous damper and a particle 
damper. The particle damper consists in a prismatic enclosure of variable 
height that contains spherical grains (total mass m p ). Contrary to what it 
has been discussed in previous experimental and simulation studies, we show 
that, for small containers, the system does not approach the fully detuned 
mass limit in a monotonous way. Rather, the system increases its effective 
mass up and above M + m p before reaching this expected limiting value 
(which is associated with the immobilization of the particles due to a very 
restrictive container). Moreover, we show that a similar effect appears in 
the tall container limit where the system reaches effective masses below the 
expected asymptotic value M. We present a discussion on the origin of these 
overshoot responses and the consequences for industrial applications. 

Keywords: Particle dampers, Granular materials, Effective mass 
1. Introduction 

Most mechanical systems, such as rotating machinery and aeronautic or 
aerospace structures, achieve damping through viscoelastic materials or vis- 
cous fluids. In general, viscoelastic materials and viscous fluids are very 
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effective at moderate temperatures (less than 250°C), but the performance 
of these is poor at low and high temperatures. Moreover, these materials 
degrade over time and lose their effectiveness. 

In recent years, particle dampers (PD) have been studied extensively for 
use in harsh environments where other types of damping are not efficient. A 
PD is an element that increases the structural damping by inserting dissi- 
pative particles in a box attached to the primary system or by embedding 
grains within holes in a vibrating structure jlj]. The grains absorb the ki- 
netic energy of the primary system and convert it into heat through inelastic 
collisions and friction between the particles and between the particles and 
the walls of the box or hole. This results in a highly nonlinear mechanical 
system. PD are effective over a wide frequency range pj . Moreover, PD are 
durable, inexpensive, easy to maintain and have great potential for vibration 
and noise suppression for many applications (see e.g. Q and j^j). 

Parameters such as the size and shape of the particles, density, coefficient 
of restitution, size and shape of the enclosure, and the type of excitation of 
the primary system, among many other features, are important in damping 
performance [4J. Thus, appropriate treatment of the PD in a given structure 
requires careful analysis and design. 

PD use a large number of small grains; therefore, its behavior is directly 
related to cooperative movements of those inside the cavity. The theoretical 
models derived from single particle systems j^] are not applicable to predict 
the performance of multi-particle systems. For more than 15 years, particle 
dynamics simulations have been used as a powerful tool for investigating the 
behavior of these types of granular systems 0, 0, 0, 0]- 

In previous works, particle dampers composed of containers of various 
sizes have been considered. In all of these works, the resonant frequency of 
the Single-Degree-of-Freedom (SDoF) system falls with respect to the un- 
damped system. This is generally attributed to the addition of the mass of 
the particles, m p . At very low excitation amplitudes, the system behaves as 
if the entire mass of the particles was attached to the primary mass, M, of 
the container (i.e. M + m p ). If the excitation level is increased, the reso- 
nant frequency gradually increases (and the damping performance increases). 
Eventually, the resonant frequency tends to the resonant frequency of the un- 
damped system. This overall behavior has been discussed in various papers 
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Yang [rjj has studied, experimentally, particle dampers under different 



conditions of excitation, different frequencies and variable gap size. The gap 
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size is the free space left between the granular bed and the enclosure ceiling 
when the system is at rest. He has found that, under some conditions, the 
system may display effective masses above M + m p or below M. However, a 
careful analysis of this phenomenon has not been carried out yet. 

In this paper, we discuss results, obtained through simulations via a Dis- 
crete Elements Method (DEM), on the resonant frequency shift of SDoF 
mechanical systems with granular damping. We show that, contrary to what 
has been discussed previously, for small containers, the system does not ap- 
proach the fully detuned mass limit in a monotonous way. Rather, the system 
increases its effective mass up and above M+m p before reaching the expected 
limiting value. Moreover, there is a similar effect in the tall enclosure limit, 
where the system reaches effective masses below the expected asymptotic 
value M. 



2. Discrete Elements Method 



In order to simulate the motion of the particles in the enclosure of a 
PD we use a DEM. This scheme, first used by Cundall and Strack lij], is 



widely used for numerical simulations of granular media [13?] . We implement 



a velocity Verlet algorithm [14| to update the positions (orientations) and 
velocities of the particles. Orientations are represented through quaternions 
to prevent numerical singularities jljj ]. 

We consider spherical soft particles. If Ri and Rj are the radii of two 
colliding particles, a = Ri + Rj — d{j is the normal displacement or virtual 
overlap between the spheres, where c?y is the distance between the centers 
of the two spheres. Under these conditions, the interaction force F p in the 



normal direction is based on the Hertz-Kuwabara-Kono model [16|, [17 • 

F n = -k n a 3/2 - JnVnV® (1) 

where k n = — v 2 )~ l is the normal stiffness (with E the Young's 

modulus, v the Poisson's ratio and = R^+RJ 1 ), 7 n the normal damping 
coefficient of the viscoelastic contact, and v n the relative normal velocity. 
On the other hand, the tangential force F s is based on Coulomb's law 



of friction [16|, |18|]. We used a simplified model in which the friction force 
takes the minimum value between the shear damping force and the dynamic 
friction. 

F s = - min (|7 s v s -v/a| , \^dF n \) sgn (v s ) (2) 
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where 7 S is the shear damping coefficient, v s the relative tangential velocity 
between the two spheres in contact and fi^ the dynamic friction coefficient. 
The sign function indicates that the friction force always opposes the direc- 
tion of the relative tangential velocity. 

The particles are enclosed in a prismatic container (the box or enclosure) 
built up of six flat walls with the same material properties as the particles 
defined through the parameters in Eqs. (pQ) and §2§. 

3. The SDoF model 

Figure [1] shows the model of our SDoF system with PD which is assumed 
to move only along the direction of the vertical z-axis. The primary system 
consists of a mass M = 2.37 kg, a spring K = 21500 Nm" 1 and a viscous 
damper with damping constant C. We have used two different values for 
viscous damping, C = 7.6 and 26.3 Nsm -1 . The undamped (C = and no 
PD) natural frequency of the primary system is fo = 15.16 Hz. 




Figure 1: (a) Model of the SDoF system with a particle damper, (b) Snapshots of the 
particles in the enclosure during a typical simulation. 

The PD is modeled as N = 250 spherical grains in a prismatic enclosure 
of lateral side L x = L y = 0.03675 m and different heights L z . The material 
properties of the particles (and walls) and the simulation parameters are 
listed in Table [TJ The gravitational field g = 9.8 ms -2 is considered in the 
negative vertical direction. Although the SDoF system can only move in the 
vertical direction, the particles move freely inside the enclosure. 
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The system is excited by the harmonic displacement of the base to which 
the spring and viscous damper are attached (see Fig. [Q. Let u(t) and z(t) 
be the displacement of the base and the primary mass, respectively. Then, 
the equation of motion for the system is given by 

Mz(t)+Cz(t)+Kz(t) = Cu(t) + Ku(t)+F pait (t), u{t) = Ucos(ut), (3) 

where F part (t) is the ^-component of the force resulting from all the interac- 
tions (normal and tangential) of the enclosure walls with the particles. The 
amplitude, U, and the angular frequency, u, of the harmonic vibrating base, 
are control parameters. 



Property 


Value 


Young's modulus E 


2.03 x 10 11 Nm- 2 


Density 


8030 kgm~ 3 


Poisson's ratio v 


0.28 


Friction coefficient /id 


0.3 


Normal damping coefficient 7 n 


3.660 x 10 3 kgs^nr 1 / 2 


Shear damping coefficient 7 S 


1.098 x 10 4 kgs^m- 1 / 2 


Excitation amplitude U 


0.0045 m 


Time step 5t 


8.75 x 10" 8 s 


Time of simulation 


13.12 s 


Particle radius 


0.003 m 


Total particle mass m p 


0.227 kg 



Table 1: Material properties of the particles and simulation parameters. 



We have obtained the frequency response function (FRF) for different 
enclosure heights L z . The initial condition for each simulation consists of a 
simple deposition of the particles inside the enclosure, starting from a dilute 
random arrangement, before applying the base excitation. 

4. Data analysis 

As shown in Table [U we have simulated the vibration of the system for 
13.12 s. After an initial transient, the system reaches a steady state. This 
steady state can display either regular or chaotic behavior, depending on the 
excitation frequency, box height, etc. In all cases, the final 10% of the time of 
the simulations has been used for the data analysis, which has proved to be 
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Figure 2: Examples of the FRF of the SDoF system with PD for different L z . Green 
squares: L z = 0.057 m, blue circles: L z = 0.1225 m and red diamonds: L z = 0.372 m. 
These results correspond to the simulations with C — 7.6 Nsm" 1 . The solid lines are fits 
of the FRF with an equivalent mass-spring-dashpot model without PD [see Eq. flljj]. 

sufficient to ensure that the steady state has been reached. We have studied 
frequency excitations in the range (0.5 — 30.0 Hz). 

We carry out a simple evaluation of the effective mass and effective damp- 
ing of the PD by fitting the FRF to a SDoF system with no particles in the 
enclosure. Other approximate methods such as the power flow method used 
by Yang present numerical instabilities and can yield negative effective 
masses (see e.g. (l9l|). 

The amplitude of the response X of a system with no PD is given by 



We carry out a least-squares curve fitting of the DEM data with Eq. (J3J). The 
values of K and U are fixed to the corresponding values in our simulations 
and C e g and M e g are fitting parameters. 

It is important to mention that for certain enclosure heights, the FRF is 
not a smooth function as the one corresponding to Eq. (00). However, the 
overall shape is well described by this fit. For small and large enclosures, 
the shape of the FRF of the PD is rather smooth and the fits are warranted. 
Examples of the quality of the fits are shown in Fig. [2] for three different 
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Figure 3: Frequency response function for the SDoF system with PD. Each curve corre- 
sponds to a different box height L z . Red squares: L z — 0.057 m, green circles: L z = 0.1225 
m and blue diamonds: L z — 0.372 m. The black continuous line corresponds to the SDoF 
system without particles in the enclosure. The black dashed line corresponds to an equiv- 
alent system with an added mass m p equal to the mass of the particles. These results 
correspond to the simulations with C = 7.6 Nsm -1 . The green circles correspond to value 
of L z for which the maximum effective damping is obtained. 

L z . An improvement of the fits for the FRF at intermediate L z for which 
some fluctuations are present could be achieved by using models with extra 
degrees of freedoms (e.g. a tuned mass damper). However, the characteristic 
frequencies of the PD are not well separated and the equivalent mass-spring- 
dashpot model constitute a good first order approximation. 

5. Results 

Examples of the FRF of the PD can been seen in Fig. [3] for a few box 
sizes L z . It is important to note that, in general, the gap size and not L z 
controls the dynamics. In this work all simulations are carried out with the 
same number of particles and particle size. Hence, the gap size is directly 
obtained as L z — h; with h = 0.039 m the approximate height of the granular 
bed at rest. For a study on the dependence of the dynamics on the number 
of particles and particle size see Ref. 0,01. 



7 



The general trends observed in Fig. [3] are consistent with previous ex- 
perimental and simulation works (see e.g. 0, jsj and [loj]). For small gaps 
(L z < 0.087 m) the FRF tends to the response of a system without particles 
but with an added mass equivalent to the mass m p of the particles. For 
very tall boxes (L z > 0.222 m), the response follows the one expected for an 
empty container. For intermediate values of L z , significant granular damp- 
ing is obtained with damped resonant frequencies intermediate between the 
two asymptotic cases. However, these intermediate values of L z yield more 
complex FRFs with the presence of more than one peak [lOj. 

5.1. Effective damping and effective mass 

A simple evaluation of the effective damping yield by the particles can be 
done by fitting the FRF to a SDoF system including only a viscous damper as 
discussed in Section HJ The effective damping, C e g is plotted as a function of 
L z in Fig. H] A clear optimum value of L z is predicted as in several previous 
works [7|, |20|. Notice that an increase of the viscous damping C leads to a 
less important influence of the PD on the effective damping C e $. This is due 
to the fact that the damping due to the dashpot of the SDoF system reduces 
the transfer of energy to the particles inside the enclosure. 

The two limiting cases (small L z and large L z ) have been explained 



111 . |21[ in terms of the proportion of the time that particles spend in 
contact with the enclosure. For small L z , particles are essentially fixed in 
their positions since the constraint imposed by the walls impede their rel- 
ative motion. Therefore, the particles behave as a simple mass added to 
the primary mass M of the system and provide no extra damping (in this 
limit the effective damping equals the viscous damping). Conversely, if the 
enclosure leaves sufficient room for the motion of the particles, the granular 
bed will reach a gas-like state (at the large excitation levels achieved near 
resonance) in which most of the time particles are in the air and only oc- 
casionally collide against the floor and ceiling of the box. In this case, the 
particles will barely influence the motion of the primary system and C e s falls 
back to the baseline, C, imposed by the viscous damper. 

In Fig. Owe plot the effective mass, M e g, obtained from the fits. The 
expected limiting masses M and M + m p are shown as horizontal lines for 
reference. Our data surveys a larger number of box heights in comparison 
with previous studies. We can see that, in contrast with the suggestion of 
previous studies, the two limiting cases are not approached in a monotonous 
way. The system reaches effective masses above M + m p as L z is decreased 
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Figure 4: Effective damping C e gf relative to the viscous damping C as a function of L^. 
The blue line corresponds to the SDoF system with C = 7.6 Nsm -1 and the green line 
to the system with C = 26.365 Nsm" 1 . The error bars correspond to the asymptotic 
standard error of the C e s best fit. 
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Figure 5: Effective mass as a function of L z . The blue line corresponds to the SDoF 
system with C = 7.6 Nsm -1 and the green line to the system with C — 26.365 Nsm -1 . 
The error bars correspond to the asymptotic standard error of the M e g best fit. 



and finally falls towards the limit value. Similarly, for large L z , an increase 
of the box height leads to effective masses below M before the limit value is 
approached. It is worth mentioning that a similar behavior has been observed 



in some experiments The author reports that the resonant frequency 
(which is simply related to our effective mass by M e g- = k/(2irf ) 2 ) can reach 
values above the expected for an empty enclosure. 

The presence of these overshoots when approaching the two limiting cases 
(zero gap size and infinite gap size) has been overlooked or received little 
attention in the past and is somewhat against intuition. In general, studies 
of PD are focused on the region of gap sizes where the maximum effective 
damping is observed. There, the effective mass is, as expected, intermediate 
between M and M + m p . However, design constraints may require a PD 
to work at gap sizes off the optimum damping and in the region of effective 
mass overshoot. In what follows, we discuss the origin of such behavior in 
the response of the PD by considering the internal motion of the granular 
bed. 
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5.2. Internal motion of the granular bed 

In Fig. [6] we plot the trajectory of the floor and ceiling of the enclosure 
over a few periods of excitation in the steady state regime. We have cho- 
sen frequencies close to the resonant frequency for each L z considered. This 
is because the effective values of mass and damping obtained by fitting are 
determined to a large extent by the frequencies with larger amplitude of mo- 
tion. In Fig. [6l we have indicated the position of the granular bed inside 
the enclosure by a band limited by the ^-coordinates of the uppermost and 
lowermost particle at any given time. Notice that such representation lets us 
have an indication of the density of the granular bed and the approximate 
time of impact with the box top and bottom walls. However, if the granular 
bed is somewhat dilute, the time of impact of the uppermost or lowermost 
particle does not coincide with the time at which the most substantial mo- 
mentum exchange happens between the grains and the enclosure. For this 
reason we also plot in Fig. [H] the total force exerted by the grains on the 
enclosure in the ^-direction (i.e. F part , see Eq. (j3J)) which not only provides 
a more precise assessment of the time of impact but of the intensity of such 
impacts. 

We recall here that the effective mass of an equivalent mass-spring- 
dashpot model is associated with the response of the PD in phase with the 
spring force, whereas the effective damping will be associated with the re- 
sponse in phase with the viscous force. 

As we mentioned, for very small gaps, the response of the system is similar 
to the response of an equivalent SDoF system where the total mass of the 
particles is simply added to the primary mass (i.e., M e g = M + m p ). This 
is due to the fact that the particles are not able to move in such reduced 
enclosures. As we slightly increase L z , particles behave as a dense pack that 
travels between the floor and ceiling of the enclosure (see Fig. E^a)). However, 
within an oscillation, the granular bed is in full contact with the ceiling or 
the floor during important portions of the time. During such periods, the 
granular bed can be considered essentially as an added mass m p . Notice 
that after leaving the ceiling (floor) the grains hit the floor (ceiling) before 
the primary system has reached its maximum displacement. The particles 
transfer momentum to the enclosure against the direction of the spring force. 
As a consequence, the effective inertia of the system during the impact is 
equivalent to a sudden mass increase. This effective mass increase exceeds 
the loss due to the short periods in which the granular bed is detached from 
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the floor (ceiling) (see the zero force segments in Fig. |5{a)). The overall 
result is an effective mass M e ^ above M + m p . 

If we further increase L z , the period of detachment within an oscillation 
increases, but the transfer of momentum at impact also increases leading 
to an overall increase of M e g (see Fig. E^b)). Eventually, the transfer of 
momentum at impact is exactly balanced by the loss of added mass due to 
the detachment periods. This makes the system render again M e g- = M + m p 
(see Fig.HD^c)). Interestingly, this crossover occurs when the granular bed hits 
the floor (ceiling) at the point of maximum displacement. 

There exist a range of values of L z for which M < M e g- < M + m p . In 
such cases the effective damping is rather high. This is mainly due to the 
fact that the granular bed hits the enclosure out of phase. In particular, 
the grains hit the base when both the primary mass is moving upward and 
the spring force is pulling upward (see Fig. Eld)). This results in a strong 
reduction of the maximum displacement of the system and in an effective 
added mass with respect to M. However, since the periods of detachment 
from the enclosure are significantly long, the average added mass due to the 
impacts is smaller than m p . 

If the impacts happens to be always at the time where no spring force 
is applied to the primary mass, then the transfer of momentum will not re- 
sult in an effective added mass. This happens when the primary mass pass 
through the equilibrium point of the spring (i.e., zero displacement). Indeed, 
we find that there exists a particular value of L z at which the granular bed 
hits the floor (ceiling) at this point and the effective mass obtained by fitting 
corresponds to M (i.e., no added mass, see Fig. EJ^e)). It is important to real- 
ize that this is also the value of L z at which the maximum effective damping 
is obtained (see Fig. land Ref.Q 0,0). Therefore, at the optimum L z 
where maximum damping is achieved, the effective mass coincides with the 
primary mass. This means that the addition of particles to create the PD do 
not affect the resonant frequency of the system if the optimum L z is chosen, 
which implies that under such conditions there is no need for compensation 
of the mass of the particles during design. 

As expected, an increase of L z beyond the optimum damping leads the 
granular bed to hit the enclosure in phase with the spring force (see Fig.EJT)). 
That is, the grains hit the floor when the spring pulls the system downward. 
The effect in the apparent inertial response is as if the system suffered a 
sudden mass decrease. Therefore, the effective mass is smaller than the 
primary mass (notice that here the granular bed is only in contact with the 
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Figure 6: Displacement and force of particles against the enclosure for different L z with 
C = 7.6 Nsm -1 . The solid lines show the position of floor and ceiling of the enclosure 
and the colored area indicates the limits of the granular bed defined as the position of the 
uppermost and lowermost particle, (a) L z = 0.042 m and / = 14.5 Hz. (b) L z = 0.057 
m and / = 14.5 Hz. (c) L z = 0.087 m and / = 14.5 Hz. (d) L z = 0.1095 m and / = 15.0 
Hz. (e) L z = 0.1225 m and / = 15.0 Hz. (f) L z = 0.147 m and / = 15.5 Hz. (g) 
L z = 0.372 m and / = 15.0 Hz. (h) Same as in Fig. [5] with arrows indicating the values 
of L z corresponding to each panel (a)-(g). 
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enclosure during the impacts). 

Much larger gap sizes lead the granular bed to expand significantly. The 
granular sample enters a gas-like state with only a few particles colliding 
with the enclosure in each oscillation. This transfers little momentum to 
the primary mass and the system presents an M e g which is close to M (see 
Fig. 

6. Conclusions 

We have studied a PD by means of simulations via a DEM. We have 
considered the effective mass and effective damping of the entire SDoF system 
by fitting the FRF to a simple mass-spring-dashpot system. In particular, 
we study the effect of the height L z of the enclosure. 

We have observed that the effective mass of the system reaches the two 
limits described in the literature for small and large enclosures. However, 
those limits are not approached in a monotonous way and clear overshoots 
appear. For small gap sizes, the system presents effective masses above the 
direct sum of the primary mass M and the particle mass m p . For large 
enclosures, the effective masses fall below M. 

We have observed that such behavior can be explained by considering 
both the period of time over which the granular bed is in full contact with 
the enclosure and the inertial effects due to the grains hitting the floor or 
ceiling in or out of phase with the spring force. 

Interestingly, we found that the value of L z at which M eS crosses M 
coincides with the optimum damping value described in the literature. Since 
the optimum L z should be simple to interpolate from the intersection with 
the horizontal M level in a plot of M e g vs L z , we suggest that such estimation 
can be a more suitable approach than the search for a maximum in a plot of 
C c fi versus L z . 

The overshoot effects described are present outside the range of maximal 
damping performance and might be considered of secondary interest in in- 
dustrial applications at first sight. However, design constrains may require 
a PD to work off the optimum damping and in one of the overshoot regions. 
In particular, for enclosures somewhat taller than the one corresponding to 
the optimum damping, one achieves effective masses only slightly below the 
primary mass. This implies that the resonant frequency is almost unaltered 
upon addition of the particles. Moreover, such values of L z achieve a re- 
markable damping (although not maximal) while still presenting a FRF with 
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a shape very similar to the one observed for a simple SDoF mass-spring- 
dashpot system. This may simplify the prediction of the behavior of the 
PD under such conditions. On the other hand, if the primary system has 
more degrees of freedom, the off-design resonant frequencies may fall in either 
overshoot regime and the side effects should be taken into consideration. 

Although we have studied a PD driven in the direction of gravity, similar 
results are expected if a horizontal setup is considered. It has been shown 
that near the resonant frequency the response of an impact damper does not 
depend on the relative direction between the motion of the system and the 
gravity Since the effective mass is largely determined by the response 
near the resonant frequency, the same general trends should be found in 
horizontally driven PD. 
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